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Abstract 

The energy-influx / entropy-influx relation in the Green-Naghdi Type III 
theory of heat conduction is examined within a thermodynamical framework 
a la Miiller-Liu, where that relation is not specified a priori irrespectively of 
the constitutive class under attention. It is shown that the classical assump- 
tion, i.e., that the entropy influx and the energy influx are proportional 
via the absolute temperature, holds true if heat conduction is, in a sense 
that is made precise, isotropic. In addition, it is proven that the standard 
assumption does not hold in case of transversely isotropic conduction. 

Keywords: energy influx, entropy influx, Miiller-Liu entropy principle, 
coldness, thermal displacement 



1 Introduction 

A key assumption in the standard continuum theory of heat propagation 
is that, for q, h and r,s the influx vectors and the external sources of, 
respectively, energy and entropy, the energy inflow (q, r) is proportional 
to the entropy inflow (h,s), the proportionality factor being the absolute 
temperature 

(g,r) = 0(M), *?>0. (1) 
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Whether or not this assumption is generally tenable has been the subject 
of considerable debate (see e.g. the papers by Miiller [8] and Liu [5j [6]). 
The issue can be taken up in the framework of one or another continuum 
theory of heat conduction: Podio-Guidugli [10] uses the framework of the 
standard theory, which leads to a parabolic heat equation, Bargmann and 
Steinmann pQ that of the so-called Type III theory of Green and Naghdi [3], 
which allows for propagation of heat waves with finite speed (we summarize 
these two theories in Section 2). In this paper, we use the latter framework, 
in a manner that differs from Bargmann and Steinmann's in a number of 
points. As in [10] for the standard theory, we prove: in Section 3, that 
for thermodynamic consistency (pQ) must hold true when heat conduction is, 
in a sense that we make precise, isotropic; in Section 4, that (pQ) may be 
violated in the transversely isotropic case. In our proofs, we exploit certain 
mathematical tools that are described in Appendices A and B. 

2 Heat conduction theories 

Thermal and deformational phenomena can be considered coupled, as is 
done, e.g., in thermoelasticity. For our present purposes, it is sufficient to 
consider the simplest instance, namely, heat propagation in a rigid conduc- 
tor. Consequently, we need not distinguish between reference and current 
configurations; we also leave all spatial dependences tacit. 

The standard continuum theory of heat conduction is based on the fol- 
lowing two laws: 

- (energy balance) 

pe = —divq + pr, (2) 

where p denotes the volumetric mass density, e the internal energy per unit 
volume, and e its time rate; 

- (entropy imbalance) 

pi] > —divh + ps, (3) 

with r\ the entropy. It is customarily assumed that ([I]) holds. With this and 
the notion of Helmholtz free-energy per unit mass: 

tp = e — $77, 

the following free- energy growth inequality is arrived at: 

ip < -rf& - i?" 1 q ■ W . 



2 



When a set of state variables is chosen - say, (i?, V??, i?) - this inequality is 
used to derive restrictions on the choice of the state functions delivering -0, 
rj, and q, in the manner devised by Coleman and Noll in their classic paper 

m- 

Green & Naghdi modified the classic path, in that they took the following 
entropy law as their point of departure: 
- (entropy balance) 

p£ = PV + div/i — p s, (4) 

where £ is the internal entropy production. As to internal energy, they ac- 
cepted the balance equation (J2J), implicitly excluding any internal production 
of energy; as to entropy and energy inflows, the classic proportionality ex- 
pressed by (JTJ). Consequently, they arrived at the following version of Q in 
terms of Helmholtz free energy: 

ip + rjti + h ■ W + £ = 0. 

Moreover, in their so-called Type III theory, they assumed that the consti- 
tutive mappings delivering ip,r), and h, depend on the following list of state 
variables: 

S := {a, a, Va, Vd}, (5) 

where the thermal displacement a, a field whose precise interpretation in 
statistical mechanics is still wanted, is somehow indirectly characterized by 
postulating that its time derivative equals an empirical temperature T: 

d(t) :=T{t)E 

As anticipated in the Introduction, Miiller and Liu questioned the gen- 
eral applicability of assumption (Q]) and proposed a thermodynamic format 

4 Here is a quotation from [4]: "The temperature T (on the macroscopic scale) is gener- 
ally regarded as representing (on thermolecular scale) some 'mean' velocity magnitude or 
'mean' (kinetic energy) 1 ^ 2 . With this in mind, we introduce a scalar a = a(X,t) through 
an integral of the form 

a= [ T(X,T)dr + a , (7.3) 

J t 

where to denotes some reference time and the constant ao is the initial value of a at time 
to. In view of the above interpretation associated with T and the physical dimension of the 
quantity defined by (7.3), the variable a may justifiably be called thermal displacement 
magnitude or simply thermal displacement. Alternatively, we may regard the scalar a (on 
the macroscopic scale) as representing a 'mean' displacement magnitude on the molecular 
scale and then T — a." 
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within which to test whether or not it fits a given constitutive class. This 
format amounts to considering the energy balance ([2]) as an internal con- 
straint to be appended to the entropy imbalance ([3]) after multiplication 
by a so-called 'Lagrange multiplier' A, a positive scalar to be constitutively 
specified: 

pf] + div/i — ps — \(pe + divq — pr) > 0, A > 0. (6) 

When the constitutive mappings involved in ([6]) depend on the Green-Naghdi 
list ([5]) of state variables, one can try and use a Coleman-Noll type proce- 
dure to see whether the multiplier A may be shown to be, if not equal, at 
least proportional to coldness [8], that is, the inverse of the absolute 
temperature. 

Granted an appropriate generalization of the format we just introduced, 
the same issue has been considered by Bargmann and Steinmann in [T] for 
thermoelastic materials. We here reach more complete conclusions than 
theirs under less stringent assumptions and by a different train of reasoning. 



3 Isotropic conduction implies proportionality of 
energy and entropy influxes 

On choosing the state variables as in ([5]), inequality ([6]) reads: 

P (d a r] - \d a e) a + p (daV - Ac^e) a + p (ckjaV - Xc\7 a e) • Vd+ 
P [d^ai] — A9v«e) Vci + divh — Adivq > 0, 

where 

div/i — Adivq = (d a h — \ d a q) ■ Va + (pah — Xd&q) • Vd+ 

+ {dvah - \ c\ a q) ■ V 2 a + (dyah - \ c\aq) ■ V 2 a. 

We apply the Coleman-Noll procedure [2] to the inequality that obtains 
by combining ([7]) and ([8]), that is, we require that it be satisfied whatever 
the local continuation (q, Va, V 2 a, V 2 d) of each admissible process. This 
requirement is satisfied if and only if 

sym(<9vo/i — A d\/ a q) 
sym(5 Vc i/i - A<9 V d<?) 
doj] — Xda£ = 0, 
dvaV - A d\/a£ = 0, 



0. 
0. 



(9) 
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and, moreover, inequality © reduces to 

P {daV ~ Xd a e)a + (p (dy a r) — Xdy a e) + d<yh — Xd&q) ■ Vd + (d a h — \d a q) ■ Va > 0. 
Conditions (J9J) x 2 are expedient to prove our main result: 

For each chosen value of the independent variables a, a, let the 
energy and entropy influxes and the Lagrangian multiplier be de- 
livered by isotropic constitutive mappings q(a, d, -, •), h(a, d, •, •), 
X(a, a, ■,■), the first two vector-valued, the third delivering pos- 
itive scalars. Then, the energy influx is proportional to the en- 
tropy influx via the Lagrangian multiplier, which depends neither 
on Va nor on Vd; 

h = Xq, A = A(a,d)>0. (10) 

Our proof is achieved as follows. For simplicity, we leave the dependence on 
(a, d) tacit, and write u for the third state variable in the list ©, and v for 
the fourth. On making use of the representation (|19p obtained in Appendix 
A, we set: 

h(u, v) = h\u + I12V, hj = hj(\u\,\v\,u ■ v) (j = 1, 2), 
q{u,v) = qiu + q 2 v, qj = qj(\u\, \v\, u ■ v) (j = 1,2), 

whence 

h - Xq = (hi - Xqi)u + (h 2 - Xq 2 )v. (11) 
With this, condition © 1 reads: 

(hi - Xqijl + (dihi - Xdiqi) u <g> u + (d 3 h 2 - Xd 3 q 2 ) v ® v+ 
(d 3 hi + 9i/i 2 |m| _1 - X(d 3 qi + diq^u^ 1 )^ sym(w <8> v + v <g) u) = 0, 

and has to be satisfied for all u and v (here di (i = 1,2, 3) denotes differenti- 
ation of a function with respect to the i-th of its arguments) . An application 
of the first of the two algebraic lemmas proved in Appendix B permits us to 
conclude that: 

hi - Xqi = 0, 

di(hi-Xqi) + (diX) qi = 0, 

d 3 h 2 -Xd 3 q 2 = 0, (12) 
d 3 hi - Xd 3 qi = 0, 
dih 2 - Xdiq 2 = 0. 



5 



Quite analogously, condition ©2 yields: 



h 2 - Xq 2 = 0, 

d 2 (h 2 - Xq 2 ) + {d 2 X)q 2 = 0, 

d 3 h 1 -\d 3 q 1 = 0, (13) 
d 3 h 2 -Xd 3 q 2 = 0, 
d 2 h 2 - Xd 2 q 2 = 0. 

Now, (fT2j) and (fTBj) imply that the Lagrange multiplier A cannot depend on 
u,v, that is to say, that 

A = A(a, d); 
with this, ip^ . and ([!% ■ (fTTI) becomes 

/i(a, d, Va, Vd) = A(a, d) q(a, d, Va, Vd), (14) 

thus, (fTU^ is established. 

A direct consequence of (|14|) is that, when two material bodies are in 
ideal thermal contact, that is, by definition, when neither the energy nor the 
entropy flux suffers a jump at a point of a common interface oriented by the 
normal field n: 

{q ■ nj = [h ■ n] = 0, 
then the Lagrange multiplier is also continuous at the interface: 

[A] = 0. 

At this point, in the words of Miiller, ". . . there is a price to pay for not 
having introduced the temperature so far in rational thermodynamics with 
Lagrange multipliers" ([S), p. 168). The desideratum is a proof that 

AoctT 1 , (15) 

with -i?" 1 the coldness; hence, in particular, that A is a universal function. 
This can be achieved, as Miiller himself proposed in [8], by considering a 
situation when a material body of whatsoever constitutive nature is put in 
ideal thermal contact with an ideal gas, for which the kinetic theory permits 
to show that (|15p indeed holds, provided $ is taken proportional to the 
empirical temperature T. We leave it to the reader to decide whether or 
not such an argument is convincing. For us, in that it relies on importing 
a result from a theory tacitly regarded as more foundational in nature, it is 
for sure suggestive of assuming (fT5"|) right away, without any need to imply, 
let alone accept, a vassalage between theories. 
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4 A counterexample to proportionality: transversely 
isotropic conduction 

In this Section, we show that assumption ([I]) is not tenable in general. 
Within the framework of a standard theory of heat conduction modified by 
the introduction of a Lagrangian multiplier, this result has been announced 
in [7] by Liu, who proposed to consider the case of transversely isotropic ma- 
terials to exhibit a counterexample to proportionality of energy and entropy 
influxes; the proof he offered is faulty, although amendable as indicated in 
[lOj . We here exploit the same counterexample within the Green-Naghdi 
theory of Type III. 

For e € V any chosen unit vector, consider the class of vector- valued 
mappings over V x V of the following form: 

/(it, v;e) = tpi u+(p 2 v+(p 3 e, ipi = (pi(\u\, \v\, u v, we, v-e) (i = 1,2,3) 

(16) 

(once again, we have left the dependence on (a, a) tacit). Clearly, each 
mapping / in this class is transversely isotropic with respect to the axis e, 
in the sense that it satisfies 

Qf(u,v)=f(Qu,Qv), Vu,v£V, VQeOrth(e), 

where Orth(e) denotes the continuous group of all rotations about the span 
of e. 

On using a representation if type (I16p for both h and q: 

h(u, v) = hi u + h 2 v + h 3 e, hi = hi(\u\, \v\, u ■ v, u ■ e, v ■ e) (i = 1, 2, 3), 
q(u, v) = qi u + q 2 v + q 3 e, hi = qi(\u\, \ v\, u ■ v, u ■ e, v ■ e) (i = 1, 2, 3), 
whence 

h - Xq = (hi - Xqi)u + (h 2 - Xq 2 )v + (h 3 - Xq 3 )e. 

Now, condition ^ 1 becomes: 

(/ii - Xqi)l + (dihi - Xdiqi) u®u+ (d 3 h 2 - Xd 3 q 2 ) v ® v+ 
(94/13 - A<94<?3) e ® e+ 

(d 3 hi + dih 2 \u\~ x - X(d 3 qi + di^M -1 )) sym(it ® v + v ® u) + 
\d4h1 + dih^u^ 1 — X(d±qi + diq^u^ 1 ^ sym(w ® e + e <8> u) + 
^4^2 + ^3/13 - A(<9 4 <22 + 9 3 g 3 )^ sym(t> <g> e + e <g> u) = 0, 
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for all u and v; condition yields a completely analogous identity. By 
applying Lemma 2 in Appendix B to each of these two identities, we recover 
relations (|12p and (|13p . as well as the following list of additional relations: 

<9 4 /i 3 - A <9 4 g 3 = 0, 
d^hi — A d^qi = 0, 
dih 3 - Xdiq 3 = 0, 

94/l2 - A <94<?2 + d 3 h 3 

and 

O5/13 - A<9 5 <? 3 = 0, 
d 5 h,2 - Xd 5 q 2 = 0, 
d 2 h 3 - \d2q3 = 0, 
d 5 hi - \d 5 qi + d 3 h 3 

whence 

hi — Xqi = 0, /12 — A(72 = 0, h 3 —Xq 3 = f with / = f(a,a), X = X(a,a). 

In conclusion, an influx discrepancy in the direction of the transverse isotropy 
remains: 

h - Xq = (h 3 - Xq 3 )e = f e, 

where 

h 3 (a, d, [Va|, |Vd|,Va • Va, S e a, S e d) = 

A(a, d) 93(0, d, |Va|, |Vd|, Va • Vd, d e a;, d e a) + /(a, d). 

Appendix A. Representation of an isotropic vector 
function of two vector arguments 

Here we prove a representation formula for isotropic vector-valued mappings 
of two vector arguments. As a premiss, we recall that a scalar- valued func- 
tion / of two vectorial arguments is isotropic, i.e., such that 

/(«,«) = f(Qu,Qv), Vu,v € V, VQ € Orth, 

where V is a suitable vector space and Orth the full orthogonal group, if and 
only if 

f(u,v) = ip(\u\, \v\,u ■ v). (17) 



- Xd 3 q 3 = 0, 



- Xd 3 q 3 = 0, 
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Lemma 1 Let / be an isotropic mapping of V into itself, i.e., such that 
Qf{u,v)=f{Qu,Qv), Vu,v 6 V, VQ G Orth. (18) 

Then, / has the following representation: 

f(u, v) = (pi(\u\,\v\,u ■ v) u + (p 2 (\u\,\v\,u ■ v) v. (19) 

Proof Write (|18p for Q a point of a smooth curve through I on the manifold 
Orth: 

t^Q(t), Q(t ) = I, Q(t ) =: W G Skw, 
and differentiate at t = to, so as to get: 

W/(u, t>) =^i(u, w) + F 2 (u, v) Wv, 

Fi(u, v) := dj{u, v), F 2 {u, v) := d 2 f(u, v), Vu, v G V, MW G Skw. 
This condition is equivalent to 

W-(a®/(tt, v)-F±(u, v)a®u-f(u, v)-Fj(u, v)a®v), Vu, v, a G V, V G Skw, 
and then 

a x f(u, v) — Ff(u, v)a x u — F 2 (u, v)a x v = , Vu, t>, a G V. 

Now, take the vector product with w := u x v and use the identity (a x 
b) x c = (a • c)b — (b • c)a so as to obtain: 

(a-w)f(u, v)— (f(u, v)-w) a = [a-Fi(u, v)w) u+(a-F 2 (u, v)w)v Vtt, v, a G V. 
Since 

(a • F±(u, v)w)u = (u ® F\(u, v)w) a = {Ff(u, v)u ® to) a, 
(a • F 2 (u, v)w) v = (v <8) F 2 (u, v)w) a = [F 2 {u, v)v ® to) a, 

we have that 

(a-w)f (it, v)—(f(u, v)-w) a = (Fi(u, v)u®w+F 2 {u, v)v®w) a, Vtt, f , a G V. 

Given the arbitrariness of a, it suffices to choose a • to = to deduce that 

f(u,v)-w = 0, Vw,t;GV <^ f(u,v) = fi(u,v)u + f 2 (u,v)v. 

To conclude, one observes that this provisional form of /(it, t;) is compatible 
with (|18|) iff fi and /2 are isotropic. Thus, an application of (|17p suffices. □ 
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Appendix B. Two algebraic lemmas 

lemma Let the algebraic equality 

al + (3u ® u + jv v + 6(u <®v + v®u) = (20) 
hold for all u, v € V. Then, 

a = /3 = 7 = ( 5 = 0. 

Proof Let it, u be chosen orthogonal and unimodular. Right-multiplication 
of (f20|) by m first and then by v yields: 

au + /3« + 5v = 0, 
av + + 5u = 0, 

whence, given that u ■ v = 0, 

a + /3 = 0, q + 7 = 0, J = 0. (21) 

On the other hand, taking the trace of (|20p gives: 

3a + /3 + 7 = 0. (22) 

To conclude the proof, it is enough to solve system (|2ip i 2 - (|22p , □ 

lemma Let the algebraic equality 

al + Bu u + <g> v + S(u <g> v + v <g> u)+ 

(23) 

ipe <g> e + x{u <8> e + e <8> u) + ^(v <S>e + e£g>i>) = 

hold for all ti, w € V. Then, 

a = /3 = 7 = 5 = (/9 = x = V ; = 0. 

Proof On exploiting the arbitrariness of u and u, we choose u ■ v = 0, 
m x v = e, |u| = \v\ = 1. Right-multiplication of (|23p by w, u, and e, 
yields: 

ait + fiu + 5v + x& = 0, 
av + 7^ + 5u + = 0, 
ae + </?e + \ u + V 1 ^ = 0, 

whence 

a + /3 = 0, a + 7 = 0, a + v? = 0, J = x = ^ = 0. (24) 
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Next, one first writes (j23|) for u = v, u ■ e = 0, then applies the resulting 
tensor to u, so as to arrive at 

(a + /3 + 7)u = 0, 

with the use of (|24() 1 . In view of the arbitrariness of u, the last relation is 
tantamount to 

a + /3 + 7 = 0; 

this, together with (|24p 1 3 . implies the desired conclusion. □ 
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